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Abstract 

In this paper we present thoughtful comments on the paper 'Bessel beams and 
signal propagation' showing that the main claims of that paper are wrong. Moreover, 
we take the opportunity to show the non trivial and indeed surprising result that a 
scalar pulse (i.e., a wave train of compact support in the time domain) that is solution 
of the homogeneous wave equation ( vector {E, B) pulse that is solution of Maxwell 
equations) is such that its wave front in some cases does travel with speed greater than 
c, the speed of light . In order for a pulse to posses a front that travels with speed c, an 
additional condition must be satisfied, namely the pulse must have finite energy. When 
this condition is fulfilled the pulse still can show peaks propagating with superluminal 
(or subluminal) velocities, but now its wave front travels at speed c. These results 
are important because they explain several experimental results obtained in recent 
experiments, where superluminal velocities have been observed, without implying in 
any breakdown of the Principle of Relativity. 

In this paper we present some thoughtful comments (Ci — C4) concerning statements 
presented in the paper 'Bessel beams and signal propagation' [1] and also some non trivial 
results concerning superluminal propagation of peaks in particular electromagnetic pulses in 
nondispersive media. 

In [1] the author recalls that the experimental results presented in [2] showed that Bessel 
beams generated at microwave frequencies have a group velocity greater than the velocity 
of light c (in what follows we use units such that c = !)[]. His intention was then to show 

"'^In [3] we scrutinized the experimental results of [2]. We presented there a simple model showing that all 
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that the signal velocity, defined according to Brillouin and Sommerfeld [B&S) was also 
superluminal. We explicitly shows that the particular example used by the author of [1], 
given by the Bessel beam of his eq.(3) does not endorse his claim. Contrary to the author's 
conclusion this beam has no fronts in both space and time domains, hence cannot satisfy 
B&S defintion of a signal. Moreover, the beam given by eq.(3) of [1] travels rigidly with 
a superluminal speed. We prove then that there are two classes of general Bessel pulses 
satisfying B&S definition of signal. A solution of the HWE corresponding to class I is such 
that the group speed is always less than c whereas its front moves with speed c.0 A solution 
of the HWE of the class II travels rigidly at superluminal speed if care is not taken of the 
energy content of the pulse. We present also some necessary comments concerning solutions 
of Maxwell equations associated with Bessel beams of classes I and II. 

We start by recalling the general solution of the HWE □$ = in Minkowski spacetime 
{M,ri,D) [10-12]. In a given Lorentz reference frame [10-12] I = d/dt & sec TM, we choose 
cylindrical coordinates (p, cp, z) naturally adapted to the / reference frame, where p = (x^ + 
1/^)2 and X = pcosip and y = psimp, with {x,y,z) being the usual cartesian coordinates 
naturally adapted to I. Writting 

^t,p,^,z) = Mp)f2{v)fsit,z), (1) 

and substituting eq.(l) in the HWE we get the following equations (where u and Q are 
separation parameters) , 



/i = 0. 



dp2 r dp 

£. + -^')h-0. . (2) 



The first of eqs.(2) is Bessel's equation, the second one implies that z/ must be an integer 
and the third is a Klein-Gordon equation in two dimensional Minkowski spacetime.0 In what 
follows ( without loss of generality for the objectives of the present paper) we choose i/ = 
(and also f2 > 0). Then, we obtain as a solution of eqs.(2) a wave propagating in the 
2;-direction, i.e., 

$jo(t,p, z) = Jo{pn) eyip[-i{ujt - kz)], (3) 



particulars of the data (including the slowing of the superluminal velocity of the peak along the propagation 
direction) can be qualitatively and quantitatively understood as a scissor's like effect. Moreover in [3] 
we called the readers attention that in [4] peaks of finite aperture approximations (FAA) to particular 
acoustical Bessel pulses called X-waves (first discoverd by Lu and Greenleaf ([5,6]) have been see to travel at 
supersonic speed i.e., with velocity greater than Cs , the sound speed parameter appearing on the homogenous 
wave equation {HWE). In [4] and [7] it is also predicted the possibilty of launching FAA to superluminal 
electromagnetic X-waves, a fact that has been confirmed experimentally in the microwave region in [2] and 
in the optical region in [8]. A review concerning the different facets of 'superluminal' wave motion under 
different physical conditions can be found in [9]. 

^Of course, this is a kind of generalized reshaping phenomena which cannot endures for ever. It lasts 
until the peak of the wave catches the front. 

■^In 4-dimensional spacetime the Klein-Gordon equation possess families of luminal and superluminal 
solutions, besides subluminal solutions. See [4] and references therein. 
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where the following dispersion relation must necessarily be satisfied, 



uj-'-e = (4) 

The dispersion relation given by eq.(4) may look strange at first sight, but there are 
evidences that it can be realized in nature (see below) in some special circunstances. 

Ci. It is quite clear that the wave described by eq.(3), called in [1] a Bessel beamQ, has 
phase velocity Vph = u/k > 1. However, we point out that the statement done in [1] is false, 
namely: 'As known, in the absence of dispersion the group velocity Vgr of a Bessel pulse 
is equal to the phase one [4,5]Q since all the components at different frequencies propagate 
with the same velocity'. To prove its falsity recall that there exists a Lorentz reference frame 
[10-12] 

/' = (1 - v^g,) ^ {d/dt + Vgrd/dz) e sec TM, (5) 

which is moving with velocity Vg^ = du/dk < 1 in relation to the frame I in the 2;-direction. 
In the coordinates naturally adapted to the frame I' the frequency of the wave is u' = Q, 
which means that in the frame /' the Bessel beam is stationary. This proves our statement 
that for Bessel beam the group velocity is always less than the velocity of light c. 

C2. Now, we show how to build two different classes (I and II) of solutions of the HWE 
by appropriate linear superpositions of waves of the form given by our eq.(3). 

Class I. Suppose, following B&S [13,14 ] that a signal is defined as a pulse with a finite 
time duration at the origin 2; = where a physical device generated it. We model our 
problem as a Sommerfeld problem [15] for the HWE (with cylindrical symmetry), i.e., we 
want to find the solution of the HWE with the following conditions (called in what follows 
Sommerfeld conditions), 

$(t,p,0) = AJo{pn)[e{t) -e{t-T)]smuJot 

1 r . Xe^-^-ij 



AJoipn)—^ [ due 
2tx 



UO — UOq 



d^{t,p,z) 



dz 



2=0 



1 r - le^^^ - 1) 
AJo{pn)—^ / duj k{uj) e-^^*^ ^. (6) 

271 J Y UJ — UJq 



In eq.(6) 0(t) is the Heaviside function, A and uoq = Vt are constants, 3? means real part 
and k{ijj) is given below and for simplicity we take T = Ntq = 2ttN/ujq, with an integer. 
Now, to solve our problem it is enough to get a solution of the third of eqs.(2). We have, 

f^(t, Z) = ^3? J J"^^ |g-*-(t-T-^s.2) _ g-^^(^-^,,,2)| 

where Vgr = k{uj)/u and V is an appropriate path in the complex a;-plane. We note 
lim^^oo'^gr = 1- Putting eq.(7) into the third of eqs.(2) we see that the dispersion rela- 
tion given by eq.(4) must be satisfied. To continue we write, 

k{u) = ^/{uj + n){uj-n). (8) 
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Figure 1: Contour for integration of eq.(7) for t — z < 



There are two branch points at ti; = ±i7. The corresponding branch cuts can be taken as 
the segment (— ^l) in the real ix;-axis. Following F from positive values of above and 
close to the real axis, the root in eq.(8) acquires a phase factor e*'" = — 1 when passing from 
^uj > Q to ^uj < —Q. Then, on the real cu-axis we have. 



a result that is necessary in order to calculate the value of /s for (t — Vgrz) > 0. We are not 
going to investigate this case here, since we are interested in the behavior of /s for the case 
where [t — z < 0). In this case, we must close the contour F in the upper half plane. Since 
there are no poles inside the contour we get that 

f^{t,z) = for t-z<0. (10) 

Now, it is easy to verify the intensity of the wave which is solution of the HWE and 
satisfies the Sommerfeld conditions given by eq.(6) has a maximum for uj = ujq, i.e., the 
waves with frequency near uq have always a much greater amplitude than all others. Under 
these conditions let us write, 

_ _ ^ 
out — kz = {ujot — koz) + {t )(ti; — ujo), (11) 

VgrO 

where Vgro = {du / dk)\uj=i^g < 1 and Vpho = ujo/ko > 1. We can write an approximation for 
the function fzit^z) denoted by f^it^z) as, 

g-«^o(t-z/fpM) j [^^-iuj(t-T-z/varo) ^-Mt-z/varo)^ ^ ^]_2) 

UJ0 — /\UJ 



"'Note that in [1] the author writes = wsin0 and k — utcosO. 

^The references [4,5] in [1] are the references [8,13] in the present paper. 
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We see that /^{t, 0) is equal to f^it, 0) if we suppress in the expression for this function the 
frequencies very different from uq. Now, fsit, 0) has support on the whole temporal axis, i.e., 
in the interval —oo < t < oo, but it is taken by some authors (like, e.g., [16]) as representing 
a wave that begin gradually at t = and ends gradually at t = T. Of course, no wave of the 
kind of /s can be build by any physical device. The importance of the function fzit.z) is 
that, as emphasized by B&S [13,14] it shows that we can associate a group velocity to pulse 
peaks in general (and of Bessel beams in particular) satisfying the Sommerfeld conditons 
(eq.(6)) and that the group velocity in this case is less than the velocity of light. This means 
that after a while the hack end of the wave that is travelling at speed c(= 1) will catch the 
peak. The wave reshapes even when propagating in vacuum. 

A general subluminal Jo-Bessel beam can be written as, 

^B{t,p,z) = Upuj)r~^[T{u)]^~^' (13) 

where T{ijj) is an appropriate transfer function and is the inverse Fourier transform. 
Now, the peaks of FAA to acoustical pulses of the form given by eq.(13) (i.e., the waves 
at 2 = are not zero only in the time interval < t < T) have been seen travelling at 
subluminal speed^ in an experiment described in [4] , thus endorsing the above analysis. 
Class II. We now return to the dispersion relation given by eq.(4) and write, 

k = kcos0, Q = ksme, (14) 

where 6* is a constant called axicon angle [5,6,17]. It results that 

Lu = ±k. (15) 

We immediately verify that 

Mtop sin e)e-'^'^'-''''°''^\ (16) 

is a solution of the HWE whose beam width is proportional to l/tu sin 9, thus being frequency 
dependent. The dependency of the beam width on frequency will cause the beam to have a 
pulse response that is independent of position. Indeed, suppose that the source is driven by 
a frequency distribution B{u), i.e., we have a pulse 



dujB{uj)Jo{ujp sin 6)e 



-i{u!t—kz cos 9) 



(17) 



If Jq were not dependent on frequency the integral in eq.(17) would be simply the inverse 
Fourier transform of the source spectrum and we return to class I solutions. However, here Jo 
is dependent on frequency and also on position and consequently modifies the pulse spectrum 
in such a way to make the time response of the pulse dependent on radial position. We put an 
index X in the wave given by eq.(17) because pulses of this kind have been named X-waves 

^Of course, in this case the speed paramenter appearing in the HWE must be Cs, the sound speed in the 
medium, and the word subluminal speed used must be understood as a speed less than Cj. 
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by Lu and Greenleaf since 1992 [5,6]. Even more, taking B(uj) = Ae"""''^' {A and Oq > 
being constants), we can easily verify (c.r., pages 707 and 763 of [18]) that we can write for 
sin^ > 0, 



/oo 
ciwe-'^ol"! Jo {ujp sin e)e-'^^'-' ™^ 
-oo 



;i8a) 



Jo 



A 



- + 



A 



[p2 sin^ e+[ao + ipYY [p2 sin^ ^ + [ao - ip] ^] 



21 I 



(18b) 



AV2{ [p^ sin^ e + al - p^Y + Aalp^ +p^smH + al-p' 



[p^ sin^ e + al - p^Y + Aalp^^' 



(18c) 

where p = {t — zcos9). 

Eq.(18c) shows that this wave is a real solution of the HWE. We recall that if in eq.(18a) 
we use as integration interval < u < oo, we get only the first term in eq.(18b). In this case 
we have a complex wave that has been called the broad band X-wave in [4-6]. These waves 
and the more general ones given by eq.(18b) propagate without distortion with superluminal 
velocity given by 1/ cos 6*, but of course they cannot be produced in the physical world 
because (like the plane wave solutions of the HWE) they have infinity energy, as it is easy to 
verify. Waves that are solutions of the linear relativistic wave equations and that propagate 
in a distortion free mode, have been called UPWs (undistorted progressive waves) in [4]. 

Now, we show that a X-pulse even if it has compact support in the time domain (thus 
being of the form of a B&S signal) is such that its front propagates with superluminal 
speed. To prove our statement we look for a solution of the HWE satisfying the following 
Sommerfeld conditions^. 



^x{t,p,0) 
d<l>xit,p, z) 



[e(t + T) -0(t-r)] / dujB{uj)jQ{ujpsme)e- 



■iujt 



dz 



[e{t + T) -e{t-T)] cose J dujB{uj)Jo{ujpsme)k{uj)e- 



(19) 



z=0 



B{lo) is taken in this example as a function such that / (hjB{uj)Jt){Lopsm9)e *'^* has support in the 

— oo 

interval —oo < t < oo. 
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and k{uj) = uo. Proceeding in the same way as in the Sommerfeld problem of class I solution 
presented above we obtain as a solution of the HWE (for 2; > 0), 



$x(t,p,^) = 7;- dujB{uj)Jo{ujpsm( 
Ztx 



^^^-i^^t-z cos 9) 



J{Lj-iD)T _ -i{u-Cj)T 



l[UJ — UO] 



-00 

CO 



/ rftu5(Lu)Jo(^psin^)e-*^(*-^™'^^) for \t- zcose\ <T 

—00 

for It — z cos 01 > T 



(20) 




Figure 2: Contours for integration of eq. (20). Fi for |t — 2;cos0| > T and F2 for |t — 2; cos 6'| < 
T. 

We see that for \t — z cos9\ < T the integral in eq.(20) is not zero. Since the axicon angle 
> 0, then 1 > cos 9 < Q and it follows that the pulse is not zero for z > t and t > T, what 
means that the wave front of our pulse propagates with superluminal speed! Of course, the 
pulse is zero for z < (t — T)/ cos 9 or z > (t + T)/ cos 9. We observe that the above result is 
true even a single Bessel pulse, i.e., when B^u) = 6{u! — uq), a result that we mentioned in 
[3]. 

How to compare this finding with the famous B&S result [13,14] stating that a wave pulse 
which propagates in a dispersive medium with loss has a front propagating at maximum 
speed c? Some things are to be recalled in order to get a meaningful answer. The first is 
that B&S example refers to a propagation of a 'plane' wave truncated in time (which, of 
course, has infinite energy) satisfying the Sommerfeld conditions (analogous to eq.(6)) and 
propagating in a dispersive medium with loss. A careful analysis [19] shows that the same 
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problem in a dispersive medium with gain reveals that in this case we can find two kinds 
of solutions ( both of of infinite energy). In one of these kinds, by appropriately choosing 
the integration path in the complex a;-plane we obtain as result that the front of the wave 
may travel with superluminal speed. This situation is somewhat analogous to what happen 
with some possible mathematical solutions of the tachyonic Klein-Gordon equation in two 
dimensional Minkowski spacetime [20,21]. This equation is important because it can be 
associated with the so called telegraphist equation. 

The reason for our finding that the X-pulse propagating in a nondispersive medium, 
although of compact support in the time domain, is such that its front travel at superluminal 
speed is the following; the solution given by eq.(20) is not of compact support in the space 
domain and as such has infinite energy as can be easily verified. Only for a pulse of finite 
energy we can warranty that its front always travel with a speed that cannot be greater than 
maximum speed. Indeed, suppose we produce on the plane z = a pulse like the one given by 
eq.(20), except that it has a finite lateral circular width of radius a, i.e., it is taken as zero for 
p > a. Such a pulse is called a FA A to the pulse given by eq.(20) and as can be easily verified 
has finite energy. If such a pulse does not spread with infinite velocity during its build up, 
then after it is ready, i.e., at t = T it occupies a region of compact support in space given by 
\x\ < R, where R is the maximum linear dimension involved. Such a field configuration can 
then be taken as part of the initial conditions for a strictly hyperbolic Cauchy problem at 
t = T. For such a problem it is well known the mathematical theorem that stablishes that 
[22,23] the time evolution of the pulse must be such that it is null for |x| > R + c{t — T). In 
conclusion, it is not sufficient for a wave to be of compact support in the time domain (i.e., 
to be a pulse) to assure that the wave front of the pulse moves in a nondispersive medium 
at maximum speed c. In order for the wave front to move with velocity c it is necessary that 
the pulse possess finite energy, and in order for this condition to be satisfied the pulse must 
have compact support in the space domain after its build up. We recall here that in [4] the 
peaks of FAA to acoustical pulses given by eq.(18) (with appropriated B{u)) have been seen 
traveling with velocities c^/ cos 9, thus confirming the theory developed above. 

C3. We now examine the claim of [1] that a wave given by our eq.(17), with B{uj) = 1, 

i.e., 

00 

U{t,p,z)= j dujJoiujpsm9)e-'^^'-'''"'''^\ u = k. (21) 
—00 

is a pulse with support only in the ^-axis at points z = t/ cos 6* and with value at that points 
(5(0). The calculations presented in [1] are wrong. Before we prove our statement let us recall 
that [1] quotes Brillouin: 'a signal can be defined as a pulse of finite temporal extension, that 
is, of infinite extension in the frequecy domain'.^ The wave given by eq.(21) has an infinite 
extension in the frequency domain but it is not a pulse of finite time domain (for a fixed 
z). Indeed, as theorem 11 on page 22 in Sneddon's book [24] stablishes: a function which 
is bounded in the time domain has an infinite extension in the frequency domain, but it is 
not true that a function with an infinite frequency spectrum is necessarily bounded in the 

^This definition is due to Sommerfeld. See [13,14]. 
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time domain. A trivial example of the last statement is the case of a Gaussian pulse, whose 
Fourier transform is itself a Gaussian. In the particular case of the wave given by eq.(21) it 
is immediate to realize that the integral is nothing more than the Fourier transform of a Jq 
function, and the value of the integral is given in many books, in particular on page 523 of 
Sneddon's book [24] . We have, 

oo 

J dujJo{ujpsiYie)e-'^'^'-^'^°'^^ (22a) 



i , = for \t — z cos 9\ < p sin 9 

Vp2sin2e-(t-2cose)2 (22b) 
for \t — 2; cos 6*1 > psin6' 

Eq.(22b) shows that U{t,p,z) has support in the entire time axis provided that \t — 
2; cos ^ I < psin^. When p = 0, since U is real (as can be seen directly from eq.(22a) we must 
have that \t — zcos6\ = and the function U is singular. We see that the result expressed 
by eq.(22b) is compatible with the one given by eq.(18b) if we take the limit for gq — > 0. 

C4. Finally, we investigate the claim (done in [1] and attributed to [8]) that the wave 
function given by eq.(3) represents an electric field. This claim is a nonsequitur. Indeed, 
the scalar solutions of the HWE can be used to generated solutions of the Maxwell system 
using the i/erte potential method (see, e.g. [25,26]). In particular, superluminal solutions of 
the HWE can be used to produce superluminal solutions of Maxwell equations [4,7,9]. If 
we choose a magnetic Hertz potential 11^ = ^j^i it is a simple exercise to show that the 
transverse electric and magnetic fields do not show any dependence on Jq. Only the 
component of the electromagnetic field configuration has a Jq dependence. Explicitly we 
have from the well known formulas E — —d/dt{V x 11^) and B — V xV x fi^ that, 

Ep = 0, E^ = iuVlJi{ytp)c-'^^'-^'\ E, = 0, 
Bp = -kQJi{Qp)e''^'''-~^'\ 



B, = n 



e 



2 2 

oj^-P = 0. (23) 

With an electric Hertz potential we obtain a solution where only the E^ component has 
a Jq dependence. As such, we conclude that the electromagnetic beams observed in [2] and 
also in [8,17] are not Jq beams. A careful analysis of the solutions of Maxwell equations in 
cylindrical symmetry shows that there are not Jq solutions representing transverse electric 
fields. The existence of only one peak observerd in the experiments done in [2] must be due 
to the Ji/p term in E^. A more detailed analysis will be reported elsewhere. 

Our conclusions are as follows: (i) our results show that the main claims of [1] are wrong 
and/or misleading and leads to equivocated conclusions concerning recent experimental re- 
sults showing superluminal motion of peaks of particular electromagnetic field configurations 
in nondispersive media; (ii) we also prove a non trivial result, namely that the condition that 
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a wave is of finite time duration is not a sufficient condition for its front to propagate with 
the speed c. It is necessary in order for the front to travel with speed c that the pulse possess 
finite energy, and thus as explained above it must (after being prepared by the launching 
device) have support only in a compact space region when ready;^ (iii) only FAA to super- 
luminal solutions of the HWE (acoustical case) and to superluminal solutions of Maxwell 
equations can be produced in nature, because only waves of this kind have finite energy. 
These FAA exhibit peaks propagating with superluminal speeds even in the vacuum, but 
since their fronts propagate with speed c this kind of phenomenom does not implies in any 
danger for the Theory of Relativity. 
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